Preprint typeset in JHEP style - HYPER VERSION 



DESY 08-174 
BI-TP 2008/39 
SFB-CPP-08-97 
HEPTOOLS 08-046 



A complete reduction of one-loop tensor 
5- and 6-point integrals 



Th. Diakonidis J. Fleischer " \ J. Gluza K. Kajda ^, T. Riemann J. B. Tausk " 

Deutsches Elektronen-Synchrotron, DESY, Platanenallee 6, 15738 Zeuthen, Germany 
^ Fakultdt fur Physik, Universitdt Bielefeld, Universitdtsstr 25, 33615 Bielefeld, Germany 
Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice, Poland 



Abstract: We perform a complete analytical reduction of general one-loop Feynman integrals 
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formalism with extensive use of signed minors is developed for the cancellation of inverse Gram 
determinants. The 6-point tensor functions of rank R are expressed in terms of 5-point tensor 
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The resulting compact formulae allow both for a study of analytical properties and for efficient 
numerical programming. They are implemented in Fortran and Mathematica. 
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1. Introduction 



At the proton-proton collider LHC and the planned e+e" collider ILC, a large number of particles 
per event may be produced. The hope is to discover one or several Higgs bosons or supersymmetric 
particles, which are typically expected to be quite heavy. The interest is also directed to the study 
of known massive particles like the W and Z bosons or the top quark. Since the production rates 
are large, a proper description of the cross-sections will typically include one-loop corrections to 
n-particle reactions, where some of the final state particles may be massive. 

The Feynman integrals for reactions with up to four external particles have been systematically 
studied and evaluated in numerous studies. We just want to mention here the seminal papers [p 
and and the Fortran packages FF and LoopTools [Q], which represent the state of the art 
until now. The treatment of Feynman integrals with a higher multiplicity than four becomes quite 
involved if questions of efficiency and stability become vital, as it happens with the calculational 
problems related to high-dimensional phase space integrals over sums of thousands of Feynman 
diagrams with internal loops. 

In this article, we will concentrate on the evaluation of massive one-loop Feynman integrals 
with n external legs and some tensor structure, 

J m'l'' n;=icr 

where the denominators Cj have indices Vj and chords qj, 

Cj = {k — qj)"^ — m'j + ie. (1.2) 

We will study in the following the cases n = 5 with R < 3 and n = 6 with R < A, and we will 
conventionally assume g„ = 0. The space-time dimension is d = 4 — 2e. 

There are several strategies one might follow. One is the reduction of higher-point tensor 
integrals to tensor integrals with less external Unes and/or lower tensor rank ||5], ^ 0, |8|] ; a second 



approach is essentially numerical [M um or semi-numerical [11, 12, 13]. A third one rests on the 



unitarity cut method [14[ 15, T^, [T^. In this case, a one-loop amplitude is evaluated as a whole, by 
using Cutkosky rules, instead of computing loop integrals from each of the Feynman diagrams. It 
is impossible to give here a comprehensive survey of recent activities, and we would like to refer 
to e.g. Jlsl , 20, 21] for recent overviews on the subject. 



Here, we will advocate yet another approach and reduce the tensor integrals algebraically 
to sums over a small set of scalar two-, three- and four-point functions, which we assume to be 
known. Whether such a complete reduction is competitive with the other approaches might be 
disputed. Evidently, this depends on the specific problem under investigation. For a study of gauge 
invariance and of the ultraviolet (UV) and infrared (IR) singularity structure of a set of Feynman 
diagrams, it is evident that a complete reduction is advantageous, and it may also be quite useful 
for a tuned, analytical study of certain regions of potential numerical instabilities. 

We have chosen a strictly algebraic approach and will rely heavily on the algebra of signed 



minors which was worked out in detail by Melrose in [22]. One of the basic observations of 



Melrose was that in four dimensions all the scalar integrals can be reduced to scalar 4-point func- 



tions and simpler ones. In [23], a representation of arbitrary one-loop tensor integrals in terms 
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of scalar integrals was derived. The representation includes, however, scalar integrals with higher 
indices Vj and higher space-time dimensions d + 21. The subsequent reduction to scalar integrals 
with only the original indices and the generic space-time dimension d is possible with the use of 



integration-by-parts identities 1 24 ] and generalizations of them with dimensional shifts. The latter 



have been derived in [25], and a systematical application to one-loop integrals may be found in 



[26].^ Basically, the reduction problem has been solved this way for n-point functions. There was 
one attempt to use the Davydychev-Tarasov reduction for the description of one-loop contributions 



to the process e+e" — > Huu [|27|], and the numerical problems due to the five-point functions were 
discussed in some detail. To a large extent they root in the appearance of inverse powers of Gram 
determinants. This feature of the Davydychev-Tarasov reduction was identified as disadvantageous 
soon after its derivation, e.g. in [^8|], where a strategy for avoiding these problems was developed. 
Besides the problem of inverse powers of the Gram determinant of the corresponding Feynman 
diagram, there are additional kinematical singularities related to sub-diagrams. This will not be 
discussed here; we refer to e.g. [28, ||, ^, ^ ||, |l^, 17 1 and references therein. 



In this article, we investigate the reduction of tensor integrals with five and six external legs 
which are of immediate importance in applications at the LHC. In Section ^ we represent tensor 
integrals by scalar integrals in shifted space-time dimensions with shifted indices. Section |3| and 
Section ^ contain our main result. In Section ^ we go one step further in the reduction of five-point 
tensors compared to [^] and demonstrate how to cancel all inverse powers of the Gram determinant 
appearing in the Davydychev-Tarasov reduction. Earlier results for tensors of rank two may be 
found in [^. Section ^ contains the reduction of tensorial six-point functions to tensorial 5-point 



functions. The corresponding Gram determinant is identically zero [ |26[ , Ig, ^, and the reduction 
becomes quite compact. Some numerical results and a short discussion are given in Section ^. 
The numerics is obtained with two independent implementations, one made in Mathematica, and 
another one in Fortran. The Mathematica program hexagon . m with the reduction formulae is 



made publicly available [31], see also [32] for a short description. For numerical applications, one 
has to link the package with a program for the evaluation of scalar one- to four-point functions, 
e.g. with LoopTools [Q, 0,0], CutTools [0, 0], QCDLoop [0]. Appendices are devoted 
to some known, but necessary details on Gram determinants and the algebra of signed minors and 
to a short summary about the reduction of dimensionally shifted four- and five-point integrals. 

2. Representing tensor integrals by scalar integrals in shifted space-time dimensions 

At first we give the reduction of tensor integrals to a set of scalar integrals for arbitrary n-point 



functions. Following [g3|, assuming here the indices of propagators to be equal to one, = 1, 
one has: 

r-d ri 

r=l 

n-1 

= -E^r^^^' (2.1) 



i=l 



'We will extensively quote from article [|2^], so we introduce here the notation (I.num) for a reference to equation 
(num) there. 
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I'd n 

J r=l 

n— 1 ^ 

= E ^f'^^^.^s'^-^^^'^^ei, (2.2) 

r=l 

n— 1 n— 1 

= - E -'.^ ^ + ^ E ( + + )^ r'^ , (2.3) 

i,j,k=l i=l 
fd n 



r=l 
n-1 

- V a^' a" nP n;, A'^^^^ 

- 1i 1j Ik Ql ^ijkl ^n,ijkl 

i,j,k,l=l 
^ n—1 

- 2 E (a'^^h'j + g'^W; + g^'qfq^ + g^'qUj + g^'q^qf + g''q^q^)n^,I^n:tf 
+ \{g'"' + g^^ g"' + g^' g"^ ) - (2-4) 

where [d+] is an operator shifting the space-time dimension by two units and 

Ad+Y,stu- /■['^+]' A 1 



p,ijk— / 11 l + 5ri+5rj+Srk-\ 5rs -<5rt -5ru - ■ 

J r=l Cr 

d'^k 



(2.5) 



where = 4 + 2Z — 2e (observe that p is the number of scalar propagators of the "p-point 



function" and that equal lower and upper indices cancel, p < n). In ([2.2|-|2.4t), the coefficients 
nij,niji. and Uijki were introduced. These stand for the product of factorials of the number of 
equal indices: e.g. nun = ^^-.niju = 3\,niijj = 2\2\, n^kk = '^InijM = l! (indices i,j,k,l all 
different from each other). Of particular relevance are the following relations for the successive 
application of recurrence relations to reduce higher dimensional integrals: 



where 



'^ijk ^ij^ijkj 

nijkl = ^ij^ijk^ijkh (2-6) 



Vij — 1 -j- dij , 
^ijk = 1 + (^ifc + Sjk, 

^ijki = ^ + 5ji + 5m. (2.7) 
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In the next step the integrals in higher dimension have to be reduced to integrals in generic dimen- 
sion. Here particular attention has to be paid to ij^^^ . Reducing the tensor integrals, this term drops 
out in general [|^, ^. 



3. Pentagons 

We start with the reduction of the pentagons. This will also provide the basis for calculating the 
hexagons as we shall see. 

3.1 Scalar and vector integrals 

For the scalar 5 -point function the recursion relation (1.31) reads forn = 5 

(d-4)i ] ir = {:] h-yi:] n o.d 





With I^^^ finite for d = 4, we have in this limit 

E^h = -A- y A"] n, (3.2) 




i.e. the scalar five-point function is expressed in the limit d — > 4 in terms of scalar four-point 
functions, which are obtained by scratching in the five terms of the sum the s^^ scalar propagator, 
respectively. This was already derived in |22], see eq. (6.1) there. 

Similarly, for the tensor integral of rank R = 1 (vector) in ( pT] ) we obtain: 

4 

^5 = E 9f ^5,i, (3.3) 



1=1 



with 



T - P - A<^+] 



i.^-^)^it'-^t{V\ n, (3-4) 

s=l \ / 5 





where again in the limit d —>■ A the I^'^'^ disappears.^ These two cases are simple and lead to a 
direct reduction to scalar integrals, without the Gram determinant Os occurring anyway. In the 
following we want to reduce tensor integrals of higher rank and show, like in §], that also in 
these cases the Gram determinant can be cancelled. 



The Is^i should not be confused with quantities introduced in Equation (2.5). 
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3.2 R = 2 tensor integrals 

The tensor integral of rank 2 can be written without a g^,y-term: 



-'5 



(3.5) 



which is obtained by replacing g^^^ by 



fj, V 




(3.6) 



(assuming qi - ■ ■ 4-dimensional and independent) and further by reducing the integrals in ( |2.2| ) 
to generic dimension. This applies in the same manner also for the tensor integrals of higher rank. 
Reducing the integrals of highest dimension, recursion relation (1.30) is used. For the 5-point 
function several cases have been worked out in : (1.41), (1.42) and (1.43). For completeness we give 
in the Appendix explicitly the cases needed in the present work. 

One remark is needed concerning the integral I^'^^ . It is known [ p6| , ^ that it always cancels 
in the end. This provides a very useful check on our calculations, which we have performed in 
every particular case under consideration. Anticipating this cancellation, we will, for the ease of 
our discussion, drop terms proportional to /I'^^' wherever they appear in the following derivation. 



With this in mind we can write for ij in (2.2) with (B.2): 



[d+? 




[d+],s 



(3.7) 



and by means of (B. 3) we obtain: 

5 

I 

-^5,jj 
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0' 



^ 1 



0\ t s\ s 



s\ t s\ 



s,t=l,Sf^i,t 



5 \ / 5 



5 \ / 5 



5 \ / 5 



5 \ / 5- 



E 



1 / s \ /is 



s,t=l,s=^i,t 



\0 s 

W 5 \ ^ 



/|*.(3.8) 



Using again, we find 



3,st rst 
3 ! 



(3.9) 
(3.10) 



s,i=l 



wiiere 



q4,s 



1 ^ 1 



=1 



s: 



3,st 



St 



s,t=l 



(3.11) 



(3.12) 



and Xff and X/^* are defined in ( |A32| ). Finally, 



1 1 



00 



'0' 



= 1 {s 



s 
s 



).'--|;(::), 



^3 



(3.13) 



In this way we have cancelled the Gram determinant for the tensor of rank 2. For later reference, 
we note that, by taking into account ( |B.4| ), we can also write 

Eqq - 



1 1 



0' 
0, 



[d+],s 
4 



(3.14) 



s=l 



3.3 R = 3 tensor Integrals 

The tensor integral of rank 3 can be written as: 

i,j,k=l 



(3.15) 
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We will now rewrite this into another representation, thereby avoiding Gram determinants ()5 in 
the denominators of the new tensor coefficients Eijk, -Boofc- 



i,j,k=l k=l 



771 Ts I Q3,st jst I Q2,stu jstu 



,t=l 



s,t,u=l 



(3.16) 
(3.17) 



According to (2.3) we have with (3.6): 



-'5,«ifc — '^tj'^tjk'^^^ijk 



5 r[d+? , 
-'5,1 "I" 



k. 



-'5,7 / \ -'5,fe 



(3.18) 



By means of recursion (|BJJ), taking into account ( |3.7| ) and keeping in mind to drop ll'^^', we have: 



0' 
k 



1-^,1 jk 



,. r[d+V,s V / 5 



5 r[<i+],s 



\/5 \/5 \/5 \/5 

s=l I \ s=l,s^i,j ( \ 



Collecting the terms proportional to Og we have with I^^ij = • • • + 2 , , ^ Eqq and ( |3.14| ): 



s = l 



0\ /s 



0\ s 



\k \k \0 

W5\/5 \/5\/5J 



r[d+ls 



5 \/5 



,0, 



^ /fc r ^ ' 



5 \ / 5 



(3.20) 



i.e. we have already cancelled one Gram determinant. We multiply ( 3.19] ) by such that we 

can make use of 



;).(:■). 



s 
k, 



+ 



s\ 



0/ \k ' 

5\/5 \/5\/5 



(3.21) 
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which will give us another factor Og. Adding all contributions, we obtain 




in 



The symbol ^ in these equations denotes a sum ^ in terms proportional to /I*", ^^^^ 
terms proportional to and Yl^s=i terms proportional to I|. Concerning the symmetrization 
in ( p.22 ), we point out that the original expression (3.18) is obviously symmetric under {i ^ j), 



while this is not explicitly seen in ( |3.22[ ) anymore. Later on, however, this symmetry will become 
apparent again. 



All terms with factors of the type ^ . j can be considered, due to ( |3.6[ ), as belonging to some 



Qfj^iy term. For other terms we have to use ( 3.21 ), which yields terms with to be cancelled. These 
are explicitly given in the coefficients of /|, /|* and i.e. ( pOj ^3% p3| ). Apart from the 



terms in the last line of (3.22) and the 



term, the remaining contributions to the coefficients 
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of I| and /|*, inserting Xfj and Xfj, can be written as 

5 



1 



s=l 



k 



s\ s 



/ 5 V-' / 



s,i=l 



s\ s' 



k 



+ 



5J ^ 

^0 s\ ft 



i \ t s 



'0 s\ t s' 



+ 



j s\ Its 



5J 



Here the following "master formula" ( Equation (A. 13) of [22] ) is of great help: 



(3.23) 



S T 

s 



s 



S T 

i s 



+ 



S T 

i , 



, r = 0,l,...5, 



(3.24) 



which yields explicitly: 



s' 




+ 



s 

1 s 



s 
s 



(3.25) 



and: 



so that ( |3.23| ) reads: 

5 

^= -E 

s=l 



j 



s \ lis 



A; 



s' 



's\ /O s 

VV5V 



+ 



s' 
s, 



(3.26) 



-E 

Next we will use : 



's\ Os\ ts' 
J AOs Ais^ 



'0 s\ t s 



I3 



(3.27) 



0^ 

A;, 



i' 
s k 



+ 



5 \ / 5 



k \0 

w 5 \ / 



(3.28) 



As trivial as this relation may look, it plays the crucial role of splitting off in order to produce 

g^^'^ terms. It might also have been written as: 



0' 
k 



s' 
k i 



+ 



5 \ / 5 



0^ 



k 



(3.29) 
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but then it would not fulfill its purpose. 



The first term at the rhs. of ( |3.28D cancels a O^, while the second term enters the g^,y-tenns, 



all of which are collected in ( pjq ). The complete coefficient of /| in ( |3.16| ) is thus given by: 



q4,S 




Finally we have to investigate the last line of ( 3.22| ), being left with the factor 



before in (3.21). The master formula (3.24) then yields 




\0, 

5 \ / 5 



as 



(3.31) 



The 



in the second term of (3.31 ) cancels and the remaining factor is antisymmetric in s and 



t, i.e. this term drops out after summation over s, t. Using again ( |3.28| ) and dropping for the time 
being the contribution to g^y terms, we finally write the coefficients of Jg* and /I*" in the following 



way, taking care of the original (i ^ j) symmetry in ( [3.22| ): 



s: 



3,st 
ijk 




(3.32) 



and 



Q2,stu 




U S t 

i s t , 
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Oj\ (u s t\ /Oi\ (u s t\ 



t s' 



{i^k) + {j ^k)} . (3.33) 



At the end we can determine the g^^i, terms from the above by collecting all terms containing factors 
of the type I ) : 



OOj 



ijk=l 



q'^q^qt (3.34) 



where the square bracket means symmetrization of the included indices, 



(3.35) 



and use has been made of (3.6). Collecting all terms of type in ( 3.22 ) we have: 



Si 



OOi 



1 ^ 1 

'2 ^ / \ 2 
s=l I s\ 



s 



s 



0/ \ j s j \s \0 j 

W5\ /5 \/5\ /5. 



'0 s' 



:,t=l 



s 



s 



0/ \ j s \ s \0 1 

\ / 5\ / 5 V / 5 \ / 5. 



s 



+- 



t s\ s 



S\ \/5V /sV /5 



+ 



t s \ / s t 



1 

- y 



s\ (^^1 u s t\ 



TStU 



(3.36) 



The following relation can be proven by multiplication with ^ J , transforming it into the relation 



for an extensional of Equation (A. 8) of |22]: 
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(3.37) 



It turns out to be useful for the simplification of the coefficients of If* and If*" in ( |3.36| ). For the 
coefficient of /g*, we apply relation ( 3.37 ) with ^ = 0. The last term on the r.h.s. of (337) is 
combined with the term on the third line of (3.36) using (3.31): 



s t 
s t 



t s' 



t s 
s 



+ 



t s' 



s' 
Os, 



t s 



s' 



t s' 



s t' 
St. 



s \ s, 

V / 5 \ / 



s r 

St, 



(3.38) 



After summation over s and t, the last term on the r.h.s. will vanish. Furthermore we apply ( |3.24 ) 
taking r = 0. 

For the coefficient of /f*" in ( 3.36D we apply relation (3.37) with fi = u. Since If*" is sym- 
metric in s,t and u, we consider the sum over all permutations of any fixed set of values of s,t and 
u. We find that 



E 



permutations I ^ 
\ S 



s V 
s t, 



t s 
s 



u s 
J\ 



s t 
s t 



't s 



t s' 
u s 



0, (3.39) 



so that the two last terms on the r.h.s. of ( |3.37| ) can be dropped in this case. Thus we have: 

' 31 







OOj 



IE 



s=l 



+ 9E 



0, 



'0 s' 



J, 



s 
s 



s 
s 



s,t=l 



1 ^ 

I E 



5 L 



0, 



'0 s' 



J. 



t s' 
Os, 



s V 

s t. 



t s' 
s 



TSt 



s,t,u=l I S 
\ S , 



U S t 

Q St, 



s V 
s t, 



t s' 
s, 



TStU 
-'2 • 



(3.40) 



Collecting all contributions, our final result for the tensor of rank 3 can be written as: 



i,j,k=l 



k=l 



(3.41) 



-13- 



771 _ \ ^ Cl4,S TS , \ ^ Q3,St jst , \ ^ Q2,StU jstU 



(3.42) 



;,t=l 



and the coefficients S^^j^, Sf-l\ Sf-l^'^ are given in ( p.30| ), ( p.32| ) and ( p.33| ) and E'oofc in ( |3.40| ). 



4. Hexagons 



The 6-point function has the nice property that the tensors of rank R can be reduced to a sum of six 
5-point tensors of rank R — 1. This property has also been derived in [0]; an earlier demonstration 



of this property, however, has been given already in [26]. The simplification in this case is due to 
the fact that Og = 0, which has extensively been discussed in p6|]. Beyond that, in our approach, 
the above results for the 5-point tensors can be directly used, thus reducing the 6-point tensors of 
up to rank i? = 4 to scalar 4- and 3- and 2-point integrals. Particularly simple results are thus 
obtained for the 6-point tensors using the results of Appendix ^ and Sections 3J and 3^. What 
was missing in [|^] is exactly this simplification, which comes with the cancellation of the Gram 
determinant 0^; see Appendix A of that paper. 



4.1 Scalar and vector integrals 



According to (1.33) we write (see [ |22[ ] and also (1.55)): 



and (3.2) now reads: 




fc = 1, . . . ,6, 



(4.1) 



E 



r 



Or) 



jrs 
J4 . 



(4.2) 



Here we see already the general scheme of reducing 6-point functions to 5-point functions: In 



general, in any signed minor ( • • • ) c: a further column is scratched, resulting in a 
in the scalar functions a further propagator is scratched. 



As in (3.3) and (3.4), with the use of (1.57), we obtain: 

5 



1=1 



and 



(4.3) 
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T[d+] 



id-5)- 



0' 



0' 



6 r[d+] 



E 

r=l 



i' 
r , 



(4.4) 



6 



While in ( [3.4| ) the first part vanishes in the limit d ^ 4, here its disappearance is due to (1.61): 

5 /„\ 



E 

i=l 



0. 



(4.5) 



Indeed ( [4.5| ) will play a crucial role for the higher tensor reduction. The resulting form in ( [4.4[ ) is 
already the generic form for the higher tensors too! Therefore it appears useful to introduce the 
vector, applying further (A.15) of [||] and (1.61): 

5 



1 



'0^ 



E 

i=l 



1 ^ 

^E 



i=l 



r 

r 



r 

k i 



k = 0,...,6, 



summing over all 5 (dependent) vectors. Vr projected on these vectors reads: 



0' 

r , 



0^ 
0, 



Sir + (qf mj 



0' 

r . 



0^ 
0, 



With this definition we can write in a compact way: 

/6" = E<^^ 



4.2 R = 2 tensor integrals 



The equation (2.2) reads in this case: 



(4.6) 



(4.7) 



(4.8) 



(4.9) 
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and by using (1.59) we have: 



0, 
2 0^ 



U \ r=l,r-^i \ / 6 



(4.10) 



6 \ / 6 

We consider the limit d —^ 4 and use (1.67): 



9 



2 ^ ( i ^ ^^^^ 



6 



(4.11) 



Writing it like in (p.5|), 



we obtain by using (3.4) 



(4.12) 



1 ^ 



0] ^.V»V6 



6 



(4.13) 



to be compared with (4.4). For completeness we specify Ej^'', which we read off from (3.4) to be 



1 



(4.14) 



and finally: 



i=l r=l,r^i 



(4.15) 



We remark that due to (4.14), E-"^ = for r = z and correspondingly this will be the case for 
all higher tensors such that limitations like r ^ i could be dropped but are convenient to keep in 
numerical programs. 



4.3 R = 3 tensor integrals 



Equation (2.3) reads in this case: 

5 



E 1^1^^k-^J-^,k + 1^(9'"^" +9''<lr + d'^'q^Vt^^ , (4.16) 



i,j,k=l 



i=l 



-16- 



and with (1.60) we have: 



^ij^ijk H^ijk - - [d-6)- 



0' 

k 



k' 
i , 



k 



0' 



0' 



6 T[d+\ j_ 



0' 
0, 



6 T[d+? 







(4.17) 



6 



The first term on the r.h.s. is eliminated due to ( |4.5D and the next two terms cancel due to ( |4.11| ). 
Taking into account relation (3.7) now reads: 



TV _,..T[d+?,r 



I r 



6 T{d+],r 



(4.18) 



As a further representation of g'^'^ we have (see (1.75)): 



I r 



r=1...6 



(4.19) 



Using again (1.57) and the definition 



i,j,k=l 



(4.20) 



we obtain: 



1 ^ 



k 



' f) \ ' — ', \ ^ 

" I r=l,rjti,j \ y 



(4.21) 



From (13.101) and (|4.19D „ reads: 



z r 



6 pr 



00 > 



(4.22) 
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so that we get: 



«J=1 r=l,r^i,j «J=1 '■=1 



I r 



(4.23) 



where in the second term we can drop the Umitation r ^ i,j since it is automatically fulfilled due 

i r 



to the numerator 



, vanishing for r = i and r = j. Thus summation over i and j is possible, 



6 



using (4.19), with a result: 



or: 



with 



i,j=l r=l,rj^i,j 



r=l 



r=l 



(4.24) 



(4.25) 



Tfiu,r 
-'5 



4.4 i? = 4 tensor integrals 



(4.26) 



The tensor integral in ( 2.4 ) contains three different integrals in higher dimension, which have to be 
reduced or to be eliminated. We begin with I^tjki using (1.26). For convenience we use x instead 
of 4: 



^J^^ n,ijkl 



r=l 



[d+Y 



E 

s=l 



/ n 

Si 



ijks n,ijk / j 

r,s=l;r^s 



/ 
r 



^ijks r » ^n,ijk 



[n + 4 - (d + 2x)] 



/ 



/' 

[d+Y,r 



/ 



r 

/c 



E ( r ) E ^■^jk-'^^n-ijks 



n.ijk 



(4.27) 



, s=l;S7^r 

The last double sum in (|4.27[), assuming all indices i,j,k to be different, reads: 



/' 
i 



/ > '^jks'^njks 



/ 



n s=l;S7^i 



'ks 



k 



n s=l;sj^k 
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/ 



(4.28) 



Now adding corresponding terms in (4.27) and (4-.28), e.g. for r = i,we get: 



/ 
i 



= l;S7^j 



Jks \ i / ^'"^^^ 



(4.29) 



^ / n s=l \ / n 



due to (1.29). In case two indices are equal, e.g. i = j ^ k, we have: 



'0 V 
i 



11 s=l;s^j \ / n 





i 



^(1 + ,5,, + 6ksK 



[d^ 



iks 



s=l 



iks ' 



s=l 



(4.30) 



like (4-.29), i.e. if two indices agree, this integral occurs only once. As final result we have: 



\ / n V / n 



r=l;r^i,j,k 



/ \ [d+]{^-i),r 
Q ^ I n-l,ijk ' 



(4.31) 



where according to ( |4.29| ) and the discussion thereafter: 



(4.32) 



and [u^lred ~ [u^]^^'* without repetition, e.g. = I ^ M 



i 



Now, making use of riijf^i = VijVijki'ijkh we see that due to (4.5) the first part in (4.31) drops 



out after insertion into (2.4). The second contribution of (4.31) yields: 

^ n— 1 n— 1 

y^I^^i' Y ''ii'^ijk[ijk]l^ed1i9j'ik- 
/ \ ij,k=l 



(4.33) 



We have: 



^ij^ijk [ijk]red = [^i^] + 
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+ 



/ 



n,ik 



(4.34) 



with the help of which ( |4.33| ) reads: 

n— 1 n— 1 



(=1 i,j,k=l 



/ 



I 
kj 



(1 + '^ifc)^„,ifc QiQk 



(4.35) 



Using (1.67) we have for cZ = 4: 

- [g^^njkll^J Qj Qk+g '^riikl^ .^ q^q^ + g ^nijP^ ..' q'^q^ | , (4.36) 

and we see that this contribution is canceled by the last three terms of the type in (p^. 

The first three terms of this type are evaluated by means of (1.59) to yield: 



"-V-^nJk 



\n 



r=l;r^i 



r "-1'^ 



(4.37) 



Inserting this into (2.4), the first part yields a vanishing contribution due to (15) . The second term 
yields, again due to ( |4~^ : 



n-l 

E{ 



i 



(x-2) 



which cancels the last term in ( |2.4[ ) and the total contribution thus reads 
1 



n-l 



'0^ 
0, 



^ij^ijk Qj Qk Qi 



i,j,k,l=l 



I \ [d+]{^-i),r 



(4.38) 



^ E {a'-^h^ + ^^''^r^ + ^?^^.r<) t (o *-Jr"" ' (4-39) 
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reducing the 6-point tensor to 5 -point tensors in lower dimensions. For further reduction we put 



explicitly n = 6 and x = 4 and write (3.18) in the form 



F -4- 



J r 
k r 



I r 



I r 



6 j[d+?,r 

"I" 



.(4.40) 



With (4.19) it is now easy to see that the square bracket in (4.40) cancels out the second part in 



(4.39) and using the definition: 



i,j,k,l=l 



we obtain: 



^6,ijkl 



1 ^ 



/ 



n \ ^ \ r / 

" I r=l,r=^i,j,k \ 7 6 



'5,ijk' 



(4.41) 



(4.42) 



Again, with ( 0T7| ) and ( pl9| ) ^ ^ reads 



^ r 



6 ipr 



I r 
k r , 



J r 
k r 



Ennh + 2- 



^-EncxA + 2- 



6 rpr 



OOii 



(4.43) 



so: 



1/ Ap 



5 6 



r 



i,j,k=l 



r=l,rj^i,j,k 



i,j,k=l r=l,r^i,j,k 



I r \ I I r \ / J i 

j r I \krl \krj 



^oofc + 2- 



OOi 



(4.44) 



and with the same argument like the one used after ( 4.23 ) we obtain the final result: 

5 6 5 6 



-'6 



E I'lj^k E v'rE.^k+g^'Y.^'k E 

i,j,k=l r=l,rj^i,j,k k=l r=l,r^k 

+9''t.^j E <^5o,+5^'E9r E ^'^5o^, 



•r 

00k 



(4.45) 
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or: 

/A'-^P = ^^PIM^V^ (4.46) 

r=l 

with: 

ir''''= E QhjQkE.jk+ E g^'-'Qt^E.lk- (4.47) 

i,j,k=l;i,j,kj^r k=l,k^r 

5. Numerical results and discussion 

In order to illustrate the numerical results which can be obtained with the described approach, we 
will evaluate a representative collection of tensor coefficients. We rely on two implementations of 
the formalism, one has been established in Fortran, and the other one in the Mathematica package 
hexagon . m. 

In the following, we denote the scalar five-point function by Eq and the scalar six-point func- 
tion by Fq. The tensor decompositions of pentagons E and hexagons F read: 

4 



S'^ = ^gfi?„ (5.1) 

i=l 
4 

E'''' = ^q'^q^E^.+g'^'^Eoo, (5.2) 

4 4 

E'"''= E QiQjlkEijk + Y.9^'"'^^^'^^^^ (5.3) 

i,j,k=l i=l 
5 

F^ = ^gfFi, (5.4) 

i=l 

E'" = E ^r^Ji^.,, (5.5) 

5 5 

p^.ux^ q'^q';qlF,J^ + Y^g^^-qlF,o^, (5.6) 

i,j,k=l i=l 

p,uXp^ E Q^ljQklfE.^.i+Y^q^qfg'P^Foo^j. (5.7) 

i,j,k,l=l «,i=l 

Please observe the difference of E^, F^ and Eq, Fq in the following. The kinematics is visualized 



in Figure ^Jj. Deviating from the first sections, we have chosen here go = in order to stay close 
to common conventions of other numerical packages. 

For the evaluation of the scalar two-, three- and four-point functions, which appear after the 
complete reduction, we have implemented two numerical libraries: 

• For massive internal particles: Looptoo Is 2.2 [IP; 
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k+Qi 



k+q. 



Figure 5.1: Momenta flow used in the numerical examples for six- and five-point integrals. 



If there are also massless internal particles: QCDLoop-1 . 4 \p5\ \ 



We observed that Looptool s may become unstable in the presence of massless internal particles, 
while QCDLoop seems to be generally slower. Our Mathematica package has an implementation 
of only Looptool s. 

For completeness, we would like to mention also other publicly available Fortran packages for 
tensor functions, which we found useful for comparisons: 

• Six-point tensors with massive internal particles: none; 

• Five-point tensors with massive internal particles: Looptools Q ; 

• Five-point tensors with both massive and massless particles: none; 



• Five- and six-point tensors with only massless internal particles: golem9 5 [|37[]. 
The two independent numerical implementations have been checked in several ways: 

• By internal comparisons of the two codes, relying on the formulae presented in this article; 
With alternative, direct representations of the tensor integrals with sector decomposition ^ 



[38] and Mellin-Barnes representations [39 



• By simplifying the numerator structures algebraically and subsequent evaluation of the re- 
sulting integrals of lower rank; 



By direct comparison with other tensor integral packages [33, |^ 



Some of the comparisons were documented in [|32|]. 



^^We used a Mathematica interface to the GINAC package sect or.decompo sit ion in order to have a convenient 
way to evaluate tensor Feynman integrals. 
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We restrict ourselves to a few phase-space points, see Tables 5J_ to 53. The first configuration 
corresponds to the reaction gg ttqq, with external momenta generated by Madgraph [41, p^ ]. 
The second configuration comes from [37], while the third is a slight modification of the first one. 
The kinematical input is completed by adding the masses of internal particles. 

We begin with massive six-point tensors. For the kinematics introduced above, we determine 

I They are complex. 



the tensor components with our Fortran pacakge as shown in Tables to | 
finite numbers. Only independent components of the tensors are shown, all the remaining ones are 
obtained by permutations of indices. 

Selected tensor coefficients of five-point tensors for the case of massive internal particles are 
shown in Table ^?7| .^ The coefficients have been compared with LoopTools 2 . 2 and indeed we 
agree. For the massive six-point functions, there is no alternative package publicly available. 

In presence of massless internal particles, we face potential infrared singularities. Then, the 
loop functions are Laurent series in e, starting with a term proportional to ^, and one has to care 
about re-normalizations compared to our basic definition |1.1[ A popular measure is [35, 37]: 

, ,4_rf rfl-2e) r d'^k 



M 



r(l + e)r2(l-e) 



ITT' 



d/2 ■ 



(5.8) 



When discussing Feynman integrals with a dependence on inverse powers of e there appears a 
dependence of their constant terms on these conventions. For convenience of the reader, the tables 
are produced with a normalization as introduced in Equation 5.8, with the choice ^ = 1. 

For the case of six-point and five-point functions with only massless internal particles, we show 



only a few sample coefficients in Table and Table p.8i which are produced with our Fortran 
package. The phase space point chosen here is defined in Table |5.2[ We checked that, within 
double precision, we completely agree with corresponding numbers produced with golem9 5. 

Finally, to complete the list of relevant results, we show also sample tensor coefficients for the 
case of both massive and massless internal particles, for five-point tensors in Table 5.10| and for 
six-point tensors in Table 5.11. For this case with mixed internal masses, there is no other publicly 
released code available. 



"^Please notice that we show here five-point tensor coefficients, while in the case of six-point tensors we have shown 
tensor components. The tensor components are representation independent and should be preferred as numerical output. 
For the five-point tensors with massive internal particles, however, we have arranged for a one-to-one correspondence 
with output of LoopTools 2 . 2, so it might be interesting to have, in this case, the tensor coefficients instead. 



Pi 


0.21774554 E+03 


0.0 


0.0 


0.21774554 E+03 


P2 


0.21774554 E+03 


0.0 


0.0 


-0.21774554 E+03 


P3 


-0.20369415 E+03 


- 0.47579512 E-F02 


0.42126823 E+02 


0.84097181 E+02 


PA 


- 0.20907237 E+03 


0.55215961 E+02 


- 0.46692034 E+02 


-0.90010087 E+02 


P5 


- 0.68463308 E-i-Ol 


0.53063195 E+01 


0.29698267 E+0\ 


-0.31456871 E+01 


P6 


-0.15878244 E+02 


-0.12942769 E+02 


0.15953850 E+01 


0.90585932 E+01 


mi = 110.0, m2 = 120.0, ms = 130.0, rrii = 140.0, ms = 150.0, me = 160.0 



Table 5.1: The components of external four-momenta for the six-point numerics; all internal particles are 
massive. For five-point functions, we shrink line 2 and fix pi + P2 Pi in order to retain momentum 
conservation. 
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Pi 


0.5 


0.0 


A A 

u.u 


A < 

U.J 


P2 


0.5 


0.0 


0.0 


-0.5 


P3 


-0.19178191 


-0.12741180 


- 0.08262477 


-0.11713105 


PA 


-0.33662712 


0.06648281 


0.31893785 


0.08471424 


P5 


-0.21604814 


0.20363139 


-0.04415762 


-0.05710657 




Pa = -{Pl +P2+P3+Pi+P5 


), mi = • • • = 


me = 0.0 



Table 5.2: The external four-momenta for the six-point numerics; all internal particles are massless. This 
set of momenta comes from [^7[]. For five-point functions, we shrink line 2 and fix pi+p2 pi in order to 
retain momentum conservation. 

To summarize, we have presented in this article tensor integrals of rank i? < 3 for five-point 
functions and of rank < 4 for six-point functions. This is sufficient for the calculation of e.g. 
four fermion production at the LHC with NLO QCD corrections. 

There are further reactions of interest which will need higher-point functions and higher ranks 
of five- and six-point functions. The details of their reductions have been left for a later investiga- 
tion. 



Pi 


0.21774554 E+Ol 


0.0 


0.0 


0.21774554 E+Ol 


P2 


0.21774554 E+0\ 


0.0 


0.0 


-0.21774554 E+0\ 


P3 


- 0.20369415 E-i-01 


- 0.47579512 E-fOO 


0.42126823 E-i-00 


0.84097181 E+00 


P4 


- 0.20907237 E+0\ 


0.55215961 E+00 


- 0.46692034 E+00 


-0.90010087 E+00 


Pb 


- 0.68463308 E-01 


0.53063195 E-01 


0.29698267 E-01 


-0.31456871 E-01 


P6 


-0.15878244 E+00 


-0.12942769 E+00 


0.15953850 E-01 


0.90585932 E-01 




mi = 0.0, m2 = 


0.0, m3 = 0.0, m4 


= 1.7430, ms = 0.0, me = 0.0 



Table 5.3: The external four-momenta for the six-point numerics; one internal mass is finite. For five-point 
functions, we shrink line 2 and fix pi + P2 ^ Pi in order to retain momentum conservation. 
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- 0.223393 E- 


-18 -i 0.396728 E- 


-19 











0.192487 E- 


-17 + ] 


0.972635 E- 


-17 


1 




- 0.363320 E- 


-17-] 


0.11940 E-17 


2 




0.365514 E- 


-17 + ] 


0.106928 E- 


-17 


3 




0.239793 E- 


-16 + ] 


0.341928 E- 


-17 




V 










0.599459 E- 


-14-] 


0.114601 E- 


-14 





1 


0.323869 E- 


-15 + ] 


0.423754 E- 


-15 





2 


- 0.294252 E- 


-15-] 


0.375481 E- 


-15 





3 


- 0.255450 E- 


-14-] 


0.195640 E- 


-14 


1 


1 


- 0.164562 E- 


-14-] 


0.993796 E- 


-16 


1 


2 


0.920944 E- 


-16 + ] 


0.706487 E- 


-17 


1 


3 


0.347694 E- 


-15-] 


0.127190 E- 


-16 


2 


2 


- 0.1 63339 E- 


-14-] 


0.994148 E- 


-16 


2 


3 


- 0.341773 E- 


-15 + ] 


0.818678 E- 


-17 


3 


3 


- 0.413909 E- 


-14 + ] 


0.670676 E- 


-15 



Table 5.4: Tensor components for scalar, vector, and rank R — 2 six-point functions; kinematics defined in 



Table 5.1 and Figure 5.1 
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1/ 


A 






\ 




U 


u 


U 


- 0.227754 b- 


-11 - 


1 0.267244 E- 


-12 








1 


0.140271 E- 


-13 - 


1 0.119448 E- 


-12 








2 


- 0.201270 E- 


-13 + 


1 0.101968 E- 


-12 


U 


U 


3 


0.102976 E- 


-12 + 


1 0.624467 E- 


-12 


A 
U 


1 


1 

1 


A 1 QOQA/l p 

U. 1 O J W'f li- 


-1/ + 


1 u. i+Z'tzy li- 


1 1 

-iz 





1 


2 


-0.131028 E- 


-13- 


i 0.610343 E- 


-14 





1 


3 


- 0.543316 E- 


-13- 


i 0.158809 E- 


-13 





2 


2 


0.181352E- 


-12 + 


i 0.141686 E- 


-12 





2 


3 


0.506408 E- 


-13 + 


i 0.1 63568 E- 


-13 





3 


3 


0.600542 E- 


-12 + 


i 0.130733 E- 


-12 




1 


1 


- 0.563539 E- 


-13 + 


i 0.178403 E- 


-13 




1 


2 


0.210641 E- 


-13- 


i 0.584990 E- 


-14 




1 


3 


0.120482 E- 


-12- 


i 0.574688 E- 


-13 




2 


2 


- 0.201 182 E- 


-13 + 


i 0.620591 E- 


-14 




2 


3 


- 0.686164 E- 


-14 + 


i 0.205457 E- 


-14 




3 


3 


- 0.447329 E- 


-13 + 


i 0.193180 E- 


-13 


2 


2 


2 


0.582201 E- 


-13- 


i 0.163889 E- 


-13 


2 


2 


3 


0.119659 E- 


-12- 


i 0.570084 E- 


-13 


2 


3 


3 


0.457464 E- 


-13- 


i 0.181141 E- 


-13 


3 


3 


3 


0.557081 E- 


-12- 


i 0.374359 E- 


-12 



Table 5.5: Tensor components for a massive rank i? = 3 six-point function; kinematics defined in Table 
and Figure 5.1. 
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V 


A 


P 
















0.666615 E-09 + 


0.247562 E-09 











1 


- 0.200049 E- 


-10 + 


0.294036 E- 


-10 











2 


0.200975 E- 


-10- 


0.237333 E- 


-10 











3 


0.645477 E- 


-10- 


0.162236 E-09 








1 


1 


-0.116956 E- 


-10- 


0.516760 E- 


-10 








1 


2 


0.160357 E- 


-11 + 


0.222284 E- 


-11 








1 


3 


0.792692 E- 


-11 + 


0.729502 E- 


-11 








2 


2 


-0.111838 E- 


-10- 


0.513133 E- 


-10 








2 


3 


- 0.681086 E- 


-11 - 


0.708933 E- 


-11 








3 


3 


-0.804454 E- 


-10- 


0.801909 E- 


-10 





1 


1 


1 


0.100498 E- 


-10- 


0.151735 E- 


-13 





1 


1 


2 


- 0.348984 E- 


-11 - 


0.195436 E- 


-12 





1 


1 


3 


-0.211111 E- 


-10 + 


0.295212 E- 


-11 





1 


2 


2 


0.357455 E- 


-11 + 


0.662809 E- 


-14 





1 


2 


3 


0.121595 E- 


-11 - 


0.807388 E- 


-13 





1 


3 


3 


0.825803 E- 


-11 - 


0.142086 E- 


-11 





2 


2 


2 


-0.958961 E- 


-11 - 


0.585948 E- 


-12 





2 


2 


3 


- 0.209232 E- 


-10 + 


0.289031 E- 


-11 





2 


3 


3 


- 0.802359 E- 


-11 + 


0.994701 E- 


-12 





3 


3 


3 


-0.102576 E-09 + 


0.378476 E- 


-10 


1 


1 


1 


1 


- 0.246426 E- 


-10 + 


0.276326 E- 


-10 


1 


1 


1 


2 


0.915670 E- 


-12- 


0.660629 E- 


-12 


1 


1 


1 


3 


0.303529 E- 


-11 - 


0.287480 E- 


-11 


1 


1 


2 


2 


- 0.822697 E- 


-11 + 


0.919635 E- 


-11 


1 


1 


2 


3 


- 0.1 16294 E- 


-11 + 


0.100024 E- 


-11 


1 


1 


3 


3 


-0.146918 E- 


-10 + 


0.183799 E- 


-10 


1 


2 


2 


2 


0.908296 E- 


-12- 


0.654735 E- 


-12 




2 


2 


3 


0.109510 E- 


-11 - 


0.100875 E- 


-11 




2 


3 


3 


0.717342 E- 


-12- 


0.557293 E- 


-12 




3 


3 


3 


0.450661 E- 


-11 - 


0.485065 E- 


-11 


2 


2 


2 


2 


- 0.245154 E- 


-10 + 


L 0.274313 E- 


-10 


2 


2 


2 


3 


- 0.318500 E- 


-11 + 


0.279750 E- 


-11 


2 


2 


3 


3 


-0.146317E- 


-10 + 


0.182912 E- 


-10 


2 


3 


3 


3 


- 0.477335 E- 


-11 + 


0.477368 E- 


-11 


3 


3 


3 


3 


-0.730168 E- 


-10 + 


0.112865 E-09 



Table 5.6: Tensor components for a massive rank i? = 4 six-point function; kinematics defined in Table 



5.1 and Figure 5.1 
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Eq 


0.702503 E-14 + i 0.170006 E-14 


El 


3.56379 E-15- 15.58904 E-16 




-7.86411 E-16 + i 1.03994 E-15 


Eoo 


-8.18587 E-ll + i 1.80354 E-11 


El23 


3.51267 E-16 + i 9.64413 E-17 


Eool 


9.38702 E-12 + 12.18811 E-11 



Table 5.7: Selected tensor coefficients of five-point tensor functions with massive internal particles; kine- 
matics defined in Table 5.1. 





e" 






Fo 


- 57.8724994 - i 9248.84583 


-3167.69411 -12981.57728 


- 1003.89197 


Fs 


-867.761166 -Fi 859.212722 


273.495904 + i 483.076108 


153.767901 


F22 


83.1234074-1271.20343 


-75.7263181 -195.1508846 


- 30.2874673 


pooo 


- 185.635891 +i 1465.754753 


487.259427 -^1525.6914058 


174.2745041 




-2.64116950-14.28827971 


- 0.8480346995 - i 0.4557274228 


-0.1450625441 



Table 5.8: Tensor coefficients Fg, F3, F22 and tensor components F*"^", F^^^^ of six-point functions; aU 
internal particles are massless, kinematics of Table 5.2. 





6° 






Eo 


202.168496 + 13211.04072 


1022.10601 +1972.027061 


309.405823 + i 0.0 


E2 


264.996441 -1303.068452 


- 96.4696846 - i 149.228472 


- 47.5008979 + i 0.0 


E33 


1780.58042 + 12914.50734 


927.71650 + 1568.572069 


180.982111 +10.0 


Eqo 


9.56327810 + i 1.61648472E-13 


4.70734562E-14 + i 2.48689958E-14 


7.10542736E-15 +10.0 


E44A 


-1035.29689-1 1422.01085 


- 452.6401 12 - i 254.226520 


- 80.9228146 + i 0.0 


Eqoi 


- 0.81227772 - i 5.68434189 E-14 


- 2.04281037 E-14 - i 2.84217094 E-14 


- 7.10542736 E-15 + 10.0 



Table 5.9: Selected tensor coefficients of five-point tensor functions with massless internal particles; kinematics defined in Table 5.2. 





6° 


lA 




Eo 


- 0.289852933 E+04 + 1 0.228935552 E+03 


- 0.945038648 E+02 + 1 0.454178453 E+02 


0.7112330546 E+01 + i 0.0 




0.168344624 E+03 -1 0.181758172 E+02 


0.4242553725 E+01 - 1 0.338838829 E+01 


- 0.6442770877 E+00 + 1 0.0 




- 0.79409571852 E+01 + i 0.5445326927 E+00 


- 0.3008645503 E+00 + i 0.9457613783 E-01 


0.1027869989 E-01 +i 0.0 


^012 


0.2472148936 E+01 - i 0.127011969 E+00 


0.9699262574 E-01 - i 0.2560545796 E-01 


-0.2331885086 E-02 + i 0.0 


^2130 


0.2733228280 E+02 - i 0.519106421 E+02 


- 0.909476582 E+01 + i 0.1744459753 E-02 


0.2112313083 E-03 +10.0 



Table 5.10: Selected tensor components of five-point tensor functions with both massive and massless internal particles; kinematics defined in Table 5.3. 





6° 








0.2403558675 E+04 - i 0.2058213187 E+03 


0.7315208677 E+02 - i 0.4276718518 E+02 


- 0.7543148872 E+01 + i 0.0 


p2 


0.1112747404 E+03 - 1 0.6809282900 E+01 


0.4419243474 E+01 - 1 0.1201033663 E+01 


- 0.1044856909 E+00 + 1 0.0 


pl3 


- 0.1014018623 E+02 + i 0.1797332619 E+01 


- 0.5914958485 E-01 + i 0.3275539398 E+00 


0.7678550480 E-01 + i 0.0 


^123 


- 0.5007216712 E+00 + i 0.4194342396 E-01 


- 0.1642316924 E-01 + i 0.7789453935 E-02 


0.1225024390 E-02 + i 0.0 


^3210 


0.1263455978 E+00 - i 0.6509987460 E-02 


0.4610567958 E-02 - i 0.1506637282 E-02 


-0.1945123881 E-03 + 10.0 



Table 5.11: Selected tensor components of six-point tensor functions with both massive and massless internal particles; kinematics defined in Table 5.3. 
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A. Gram determinants and algebra of signed minors 



In this section relations are derived, which will turn out to be indispensable in our tensor reductions. 
We begin with some notational remarks on Gram determinants Gn-i, 



Gn-i = \2qjqk\, j,k = 1,- ■ ■ ,n- 1. 
The modified Cayley determinant of a diagram with n internal lines with chords qj is: 

On = \Cjk\ , j,k = (),■ 



, n, 



(A.l) 



(A.2) 



1 1 

1 Yn Yi2 
1 Yi2 Y22 



1 Yin Y2n 



1 

Yi, 

Y2r 



Y 



with 



Yjk = -{qj - qk) + ruj + m^. 
From our choice qn = 0, it follows that both determinants are related: 

On = —Gn-l, 



(A.3) 



(A.4) 



and we will usually call ()„ the Gram determinant of the Feynman integral. 

Signed minors [ p2[ ] are determinants (with a sign convention) which are obtained by excluding 
rows and columns from the modified Cayley determinant ()„. They are denoted by the symbol 



Jl J2 ■■■Jm 

ki k2 • • • kjn, 



(A.5) 



labelling the rows j'l, j2, ■ • • ,jm and columns ki,k2, - ■ ■ ,km which have been excluded from ()„. 
The sign of a signed minor is defined by 

(_l)ii+j2+- +i™+fci+te+...+fc„ ^ Signature[ji,i2, • ■ ■ jm] x Signature[A:i, ^2, • • • A:n^](A.6) 

where Signature gives the sign of permutations to place the indices in increasing order. This 
agrees e.g. with the definition of the operator Signature[List] in Mathematica. As an example 
may serve the quantity A„: 



1^11 1^12 . 
Y12 Y22 . 

Ym Y2n . 



Yin 
Y2n 

Y 




(A.7) 



We now will derive two relations between signed minors. Let us introduce 




(A.8) 
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We are going to siiow that this expression can be factorized as 



A' 



(A.9) 



and provide an explicit expression ioxXf-. To begin with, we show that A^- is symmetric in the 



indices i and j for fixed s. Obviously the third term on the right hand side of (A. 8) is symmetric 
since we consider a symmetric determinant. The symmetry of the first two terms means 



0\ /Oj 



0\ /O i 



i I \^ s \ i \0 s 
K / 5 \ 75 V/5\ / 5. 



+ 







S 



s 



I \ 1 s \1 \i s 

K / 5 \ 7 5 V / 5 \ 7 5J 



The first square bracket of (A.IO) can be evaluated using (A.13) of [22], i.e. 



0\ /Os 



+ 



;).(::). 



and (A.IO) then results in 

{;).(::). 



0. 



This is proved by multiplication with and using Eqn. (A.8) of [22] with r = 2, i.e. 
I i l\ I \ I i\ I l\ I i\ I I 



\k 

5 \ 7 5 



y 7 5 V/ 5 



; i,j,k,l — 0, . . . ,5. 



= 0. 
(A.IO) 

(A. 11) 



(A. 12) 



(A.13) 



Inserting this, products of three factors of the form I ^ 1 cancel by pairs, q.e.d. . 



For the following, relations (A.ll) and (A.12) of [22] become important, i.e 



E 

1=1 



()5 



and 



(A. 14) 



J 



0, (j ^ 0). 



(A. 15) 



i=l \ 7 5 

Further, "extensionals" are needed, i.e. relations valid for Os can be extended to any minor of Os; 



an extensional of (A.14) e.g. is 



j 



J 



y , . , , 

^ \ k i j \k 

i=l \ 7 5 \ 7 5 



(A. 16) 



'Assuming here ()g 7^ means no limitation since we are just looking for an algebraic relation. 
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As the simplest case we now immediately obtain from (AJ) A^^ = 0, i.e. 

XL = 0. 



(A. 17) 



Applying ( |A.14| ) and ( |A.15| ) to ( |A.8[ ), we see 

5 



(A. 18) 



and due to the symmetry in i and j we also have 

5 / \ 



1=1 



(A. 19) 



which gives us a hint of how Xfj might look, namely due to ( |A.18| ) it should contain a term 



i \ i s 



A further contribution must vanish after summing over i. Due to (AJ/7) it 



must contain a factor 



si 



The second factor of this contribution can only be depend on s 



s 



and has been determined by explicit calculation to be I I . Thus we conclude: 



X. 



s' 
i , 



s' 

5V V5 



+ 



J 

s i 



5 V / 5 



(A. 20) 



We come now to the second relation between signed minors. While ( IA.20D will be needed for 
the reduction of 4-point tensors to scalars /|, for the reduction of 3-point tensors to scalars /|* we 
also need 







t s 



i I \0i/ \s j 

V/5\ /5\/5 



5VV5 



vst 
ij ' 



(A.21) 



where again we have to show that indeed Og factorizes and we have to give an explicit expression 
for X^j. The left-most term on the left hand side is an auxiliary term. It is antisymmetric in s and t 



after the cancellation of 



and vanishes after summation over s and t because /|* is symmetric 



in s and t. The cancellation of has to be checked explicitly in every case where (A.21) is 

applied. 
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We observe that the expressions for X^j (A^) and X|* ( A.21 ) differ only by replacing one 



by t. Therefore the following ansatz is implied for X^j. 



St 



s 



t s' 

\i SI 

5 \ / 5 



+ 



i 

si, 



5 \ / 5 



(A.22) 



Now we directly evaluate XfjO^ using (A. 13) 



0' 



5VV5 



5VV5. 



and the remaining equation to be verified is 



+ 



0' 



5VV5 



0' 



t s 
s. 



(A.23) 



0, 



5VV5 



t s' 
i s . 



5VV5 



t s 
s 



'0' 



t s 
i , 



(A.24) 



which is done by multiplying again with Os and again using ( A.13 ). This gives us at the same time 
also a more general proof for Xf^ ( |A.20D , putting t = 0. 



B. Reduction of dimensionally shifted five- and four-point integrals 



In this appendix we provide explicitly the needed recursion relations for the reduction of the five- 
and four-point functions. In spite of the fact that here, essentially, only two different relations 



of [26| are applied for different indices and dimension, namely (1.30) and (1.31), we consider it 



helpful and sometimes even necessary, to provide them in detail. A special case of (1.31) is (B.4). 
The others are special cases of (1.30). For the six-point function relation (1.26) plays a major role 
and will be quoted when applied. 



0' 

k 



5 rld+]^ 



+ E ■ 



k 



5 r[d+?,s 



k 



5 T[d+? , 



k, 



^I^^(B.I) 



0' 
J. 



5 r[d4 



+ E ■ 



J. 



5 T[d+], 
4,i 



+ 



5 T[d+] 
-'5 • 



(B.2) 



The four-point function's shift is (1.44): 



s' 

1 s, 



4,i 



t S 

i s 



_5jst 



3 ) 



(B.3) 
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and the four-point integrals occurring in the reduction are (1.50): 



r[d+],s 




1 



d-3 



(B.4) 



In applications we can put d = 4 since I^^'^ is UV- and IR-finite. Beyond that, as it is done 



frequendy [|], +^ can be used as well as a "master integral" (see e.g. ( p.l4[ )) without reduction 
to the generic dimension. 
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